Persistence of a Particle in the Matheron-de Marsily Velocity Field 
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We show that the longitudinal position x(t) of a particle in a (d+ l)-dimensional layered random 
velocity field (the Matheron-de Marsily model) can be identified as a fractional Brownian motion 
(fBm) characterized by a variable Hurst exponent H(d) = 1 — d/4 for d < 2 and H(d) = 1/2 for 
d > 2. The fBm becomes marginal at d — 2. Moreover, using the known first-passage properties of 
fBm we prove analytically that the disorder averaged persistence (the probability of no zero crossing 
of the process x(t) upto time t) has a power law decay for large t with an exponent 9 — d/4 for 
d < 2 and 9 = 1/2 for d > 2 (with logarithmic correction at d — 2), results that were earlier derived 
by Redner based on heuristic arguments and supported by numerical simulations (S. Redner, Phys. 
Rev. E 56, 4967 (1997)). 



PACS numbers: 05.40.-a, 02.50.Ey, 05.60.-k 

The Matheron-de Marsily (MdM) model, originally 
introduced to study the hydrodynamic dispersion of a 
tracer particle in porous rocks [1], provides perhaps the 
simplest setting where anomalous super-diffusion arises 
due to long range correlations generated dynamically by 
the process itself. In the MdM model a single particle dif- 
fuses in a (1 4- l)-dimensional layered medium with one 
transverse and one longitudinal direction. The motion 
of the particle is purely Brownian along the transverse y 
direction. In contrast, along the longitudinal x direction 
the particle is driven by a drift velocity v(y) that is a 
random function of only the transverse coordinate y (see 
Fig. 1). Even though the velocities of different layers 
in the MdM model are uncorrelated, the motion along 
the x direction gets correlated in time due to the mul- 
tiple visits to the same transverse layer by the particle 
in a given time t. This generates a typical bias in the 
x direction giving rise to a super-diffusive longitudinal 
transport: y /f {x T Jt)) ~ £ 3 / 4 for large t [2-5]. 



FIG. 1. A realization of the MdM model in (1 + 1) di- 
mensions. The heavy black dot represents the particle. The 
velocity field, uniform in a row but fluctuating from row to 
row, is represented by the arrows. 

While the transport properties in the MdM model are 



by now well understood, the first-passage properties are 
much less obvious [4,6]. The first-passage property of 
a particle diffusing in a random medium is a subject of 
long standing interest [3,7]. Quite generally, the first- 
passage probability F(t) is simply the distribution of the 
time when a particle first crosses a given point in space. 
The persistence P(t), on the other hand, is the prob- 
ability that the particle does not cross a given point 
upto time t. The two probabilities are simply related: 
F(t) = —dP(t)/dt [7]. The calculation of the persistence 
P(t) is nontrivial even in pure systems without disor- 
der. There have been tremendous theoretical and exper- 
imental activities over the last decade to understand the 
persistence properties in nonequilibrium systems without 
quenched disorder [8]. Not surprisingly, the persistence 
in disordered systems is even harder to understand. The 
one dimensional Sinai model [9] , where a particle diffuses 
in a random potential that itself is the trace of a random 
walk in space, is perhaps the only model with quenched 
disorder for which the persistence properties are under- 
stood analytically [10,11]. It is, therefore, important to 
find other models with quenched disorder for which the 
persistence can be computed analytically. Given its rel- 
ative simplicity, the MdM model seems to be a likely 
candidate for this purpose. 

The persistence properties in the MdM model were first 
studied by Redner [4,6]. Perhaps the simplest question 
is: What is the probability P(t) that the particle does 
not re-cross the y axis upto time t, starting initially at 
x(t = 0) = 0? This is equivalent to restricting the par- 
ticle in the semi-infinite geometry (x > 0) with an ab- 
sorbing boundary condition at x = and P(t) then rep- 
resents the survival probability of the particle [7]. The 
probability P(t), of course, varies from one realization of 
the velocity field to another. Redner argued physically, 
supported by numerical simulations, that the disorder 
averaged survival probability P(t) decays as a power law 
for large t, P(t) ~ t -1 / 4 [6]. Extension of this argu- 
ment to a (d+ l)-dimensional generalization of the MdM 
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model (with d transverse directions and one longitudinal 
direction) suggests a rather rich asymptotic behavior [6] , 



in = Vi{t), 



(3) 




< d < 2 
d= 2 
d> 2, 



(1) 



indicating that d = 2 is a 'critical' dimension below which 
the particle survives longer (compared to the ordinary 
diffusion) due to the random velocity field. The argu- 
ments used by Redner in deriving these asymptotic re- 
sults, though physically intuitive, are heuristic. Hence, a 
rigorous derivation would be important and welcome. 

The main result of this Rapid Communication is to 
show that the stochastic process x(t), representing the 
longitudinal coordinate of the particle in the (d + 1)- 
dimcnsional MdM model, can be identified as a fractional 
Brownian motion (ffim) introduced by Mandelbrot and 
van Ness [12] and is characterized by a Hurst exponent 
H{d) that depends only on the dimension d. This result 
shows that the MdM model provides a physical realiza- 
tion of a fBm with a variable (by varying the dimension 
d in the MdM model) Hurst exponent. This connection 
to fBm also allows us to prove the asymptotic results 
in Eq. (1) by using a known result for the first-passage 
probability of fBm. 

For later purposes it is useful to summarize the defini- 
tion and the known first-passage properties of a fBm. A 
stochastic process x(t) (with zero mean E [x(t)] = 0) is 
called a fBm if its incremental two-time correlation func- 
tion C(h,t 2 ) = E [(x(ii) - x{t 2 )f 

depends only on the difference \t\ — 
grows asymptotically as a power law, 



is stationary, i.e., 
t 2 \ and moreover, 



C(t u t 2 ) ~ |*i -t 2 \ 2H , \h-t 2 \»l. 



(2) 



The parameter < H < 1 is called the Hurst exponent 
that characterizes the fBm [12] and E [. . .] denotes the 
expectation value over all realizations of the process. For 
example, the ordinary Brownian motion is a fBm with 
H = 1/2. The zero crossing properties of fBm has been 
studied before [13-16] and it is known, both analytically 
and numerically, that the probability that a fBm does not 
cross zero upto time t has a power law decay, P(t) ~ t~ e 
for large t with 6 = 1 — H [15,16]. Note that the result 
9 = 1 — H docs not require the process x(t) to be Gaus- 
sian and holds for any zero mean process that satisfies 
Eq. (2). Our strategy, in this paper, would be to first 
establish that the transverse coordinate x{t) of the parti- 
cle ina(d+ 1)- dimensional MdM model indeed satisfies 
Eq. (2) with a Hurst exponent H (d) and then, to use the 
known fBm result 9(d) = 1 — H (d) to prove the results 
in Eq. (1). 

In the (d+l)-dimensional MdM model considered here, 
the d > transverse coordinates of the particle, denoted 
by j/i(i)'s with i = 1, 2 . . . , d, undergo ordinary Brownian 
motion, 



where r)i 's are standard zero mean Gaussian white noises 
with correlators, (r}i(t)rjj(t')) = 8ij8(t — t'). The longitu- 
dinal coordinate x(t), in contrast, is driven by a random 
drift ^[y(i)] that depends only on the tranverse coordi- 
nates y(t) = {y t (t)}, 



x = v\y(t)]+£(t), 



(4) 



where £ (t) is again a delta correlated zero mean Gaussian 
white noise and is uncorrelated to the noises ?7i(i)'s. The 
velocity field v[y] is quenched, i.e., for a given realization 
of the function v, one first evolves x(t) by Eq. (4) and 
then one needs to 'disorder average' (denoted by ~) over 
different realizations of the random function v. In the 
original MdM model in (1 + 1) dimensions, the random 
function v(y) is considered to be a Gaussian with zero 
mean, v(y) = and a correlator v(yi)v(y 2 ) = S(yi — y 2 ). 
In the generalized (d+ l)-dimensional version considered 
here, we still choose v[y] to be a zero mean Gaussian 
function, but with a short-ranged correlator, 



v[yi}v[y 2 ] 



[2^da 2 ] d ' 2 



exp 



(yi -Y2) 
2da 2 



(5) 



The cut-off a in Eq. (5) represents the correlation length 
of the velocity field in the transverse directions. Physi- 
cally, this mimics the fact that the velocity layers have a 
finite thickness of length a. For d < 2, one can safely take 
the limit a — ► and recover the delta correlated disorder 
of the original MdM model. However, a nonzero cut-off 
is necessary for d > 2 in order that the MdM model in 
continuum space is well defined. The physical reason for 
this nonzero cut-off is due to the fact that in d > 2, a 
point Brownian particle never meets another fixed point 
in space. The events, where two point particles meet, 
constitute a set of measure zero. Hence, to obtain phys- 
ically meaningful results, it is necessary to have a finite 
size for either of the particles in d > 2. In the present 
context, it means that the point particle in our model 
will feel the velocity fields in the d transverse layers pro- 
vided the layers have a nonzero thickness. The choice of 
a Gaussian function in Eq. (5) just makes the subsequent 
calculations easier, but in principle one can choose any 
short-ranged function in Eq. (5). The exponent 9 turns 
out to be independent of this choice. 
Integrating Eq. (4) we get, 



'(*)= [ Z(r)dT+ fv[y(r)]dT, 
Jo Jo 



(6) 



where we have assumed x(0) = 0. In order to relate 
the process x(t) in Eq. (6) to a fBm, all we need is 

to compute the expectation value E {x(t\) — x(t 2 )) 



2 



where E denotes an average over all realizations of x(t), 
arising from the thermal noises as well as the disorder, 
i.e., £[••■] = (~) = (• • •). Clearly, from Eq. (6) one 
gets E [x{t)\ = 0. Using <^(*i)^(t 2 )> = <5(*i - t 2 ), one 
gets from Eq. (6), 



E [xit^xih)} = min(ti, t 2 ) + I(h,t 2 ), 
where I{t\,t 2 ) is given by 



(7) 



I{ti,t 2 ) 



ti rU 



(v [y (n )]v [y (r 2 )] )dri dr 2 



o ^0 

ti /-t 2 



= a / / ( e -(y(^)-y(r2)) /2da ) rfTldT2j (g) 
'o Jo 

with a = [2Trda 2 ]~ d / 2 . The Gaussian correlator in Eq. 
(5) has been used in going from the first to the second 
line in Eq. (8). By integrating Eq. (3), it follows that 
the vector 



FT 1 

Yi = yi(n) - j/i(r 2 ) =/ i]i{T)di 

J To 



(9) 



is a Gaussian random variable since 77^ 's are Gaussian. 
The vector Y has zero mean, (Y) = and a variance 
(Y 2 ) = d\r\ — T2I, computed easily from Eq. (9) using 
the property (Vi{ T i)Vj( T 2)) = SijS(ri — t 2 ). The quantity 
(e~ Y l 2da ), that appears inside the integral on the right 
hand side of Eq. (8), can then be evaluated using the 
explicit Gaussian form of the distribution 

P(Y) - Mn - r 2 |]- d / 2 exp [-Y 2 /2rf|n -r 2 |] , (10) 

and we get 

dT\dT2 



I(h,t 2 ) = [2ird]- d / 2 I' f 
Jo Jo 



\n-T 2 \ + a 2 } d / 2 ' 



(11) 



Fortunately, the double integral in Eq. (11) can be ex- 
plicitly performed. Using this explicit form in Eq. (7) we 
get, for all t x > 0, t 2 > and d ^ 2, 4, 

E [x(t{)x{t 2 )\ = Amm(t u t 2 ) + B [(h + a 2 ) + (t 2 + a 2 f 
- (|*i - t 2 \ + a 2 f - a 2 ?} , (12) 

where (3 = (4 - d)/2, A = 1 - 4a 2 - d (27rd)- d / 2 /(2 - d), 
and B = 4(27rd)- d / 2 /(2 - d)(4 - d). In the limit d -> 4, 
-B diverges but /3 — > 0, and the second term in Eq. (12) 
just becomes a logarithm, but stays finite. The other 
marginal case d = 2 will be discussed later. Note that 
while B is independent of the cut-off a, A depends ex- 
plicitly on a. 

Putting t\ = t 2 = t in Eq. (12), we get the results for 
the variance in the longitudinal coordinate of the particle 
for all t > 0, 



This result in Eq. (13) clearly emphasizes the role of the 
cut-off a and the critical dimension d = 2. For d < 2 or 
equivalent^ (3 = (4 - d)/2 > 1, it follows from Eq. (13) 
that for large t the second term on the right hand side 
dominates, giving rise to a super-diffusion E \x 2 (t)\ w 
2Bi (4-d)/2 Thc cut _ off fl playg nQ role for d < 2 an d in 

fact, one can safely take the a — > limit in Eq. (13) for 
d < 2. In contrast, for d > 2 (/3 = (4 - d)/2 < 1), thc 
first term on the right hand side of Eq. (13) dominates for 
large t and one gets ordinary diffusion, E \x 2 (t)\ w At 
where A = 1 - 4a 2 - d (27rd)- d / 2 /(2 - d) > depends ex- 
plicitly on the cut-off a. In fact, A diverges in the limit 
a — > 0. Thus, to get a physically meaningful result, it is 
necessary to have a finite nonzero cut-off a for d > 2. 
Thc incremental two-time correlation function 



C(ti,t 2 ) = E \{x{h) - x{t 2 )) 
Eq. (12) for all h,t 2 > 0, 



£ [x 2 (t)] = At + 2B [(t + a 



can be computed from 



C(h,t 2 ) = A |ti - i 2 | + 2B [(\h -t 2 \+ a 2 f - a 



261 



(14) 



Thus, in thc limit \ti — t 2 \ >> 1, the function C(t\,t 2 ) 
again has two different asymptotic behaviors depending 
on whether d < 2 or d > 2. In the former case, since 
13 = (4-d)/2 > 1, one gets, C(ii,i 2 ) « 2^1^ - t 2 f for 
large |ii — t 2 \. As in the case of the variance in Eq. (13), 
the cut-off a does not play any significant role for d < 2 
and one can take safely the limit a — > in Eq. (14). The 
incremental correlator is thus stationary and moreover, 
grows as a power-law for large \t\ — 1 2 \ as in Eq. (2). This 
proves that for d < 2, the longitudinal displacement x(t) 
of the particle is indeed a fBm with a Hurst exponent, 
H(d) = (3/2 = 1 — d/4. In the complementary case d > 2, 
the first term on the right hand side of Eq. (14) domi- 
nates for large \ti — 1 2 \ and one gets C(t\, t 2 ) w A\t\ — 1 2 \ 
with A > 0. Once again, as in Eq. (13), A depends ex- 
plicitly on the cut-off a (and diverges as a — > 0) and it 
is necessary to retain a nonzero cut-off a > for d > 2. 
With a finite cut-off, C{t\,t 2 ) again has the fBm form 
as in Eq. (2) with a Hurst exponent H{d) — 1/2 for all 
d > 2. This indicates that x(t) asymptotically becomes 
an ordinary Brownian motion (or fBm with H = 1/2) for 
all d > 2. This shows that for all d 7^ 2, the longitudi- 
nal coordinate x{t) of the particle is a fBm with a Hurst 
exponent, 



H(d) 



-ft; 



4 ■ 



d < 2 

d>2. 



(15) 



2ys _ a 2/3n 



(13) 



It then follows immediately from the known first- 
passage property of fBm mentioned earlier that the prob- 
ability, that the process x(t) does not cross zero upto time 
t, decays as a power law, P(t) ~ t~ e with 6 = 1 — H. 
Using the results for H in Eq. (15), we get = d/4 for 
< d < 2 and 8—1/2 for d > 2, thus proving the results 
in Eq. (1) for d ^ 2. 
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In the marginal case d = 2, one can again evaluate the 
integral I(t\,t 2 ) in Eq. (11) explicitly. Substitution of 
this form in Eq. (7) gives for all *i, t 2 > 0, 

E [x(ti)x(t2)\ = A' minih, t 2 ) + B' [(h + a 2 ) ln(*i + a 2 ) 
+ (t 2 + a 2 ) ln(i 2 + a 2 ) - a 2 ln(a 2 ) 
- (|*i - t 2 \ + a 2 ) In(|ti - t 2 \ + a 2 )] , (16) 

where A 1 = 1 - ln(a 2 e)/27r and B' = 1/4tt. For h = t 2 = 
t, one gets from Eq. (16) the exact result for the variance 
for all t, 

E [x 2 (t)] =A't + 2B' [(t + a 2 ) ln(t + a 2 ) - a 2 ln(a 2 )] . 

(17) 

Thus, asymptotically for large t, the longitudinal mo- 
tion for d = 2 is super-diffusive, but only marginally 
E [x 2 (t)] « 2B't\n{t). One can easily com- 
pute the incremental two-time correlator C{t\,t 2 ) = 



E 



(x(h) - x{t 2 )) from Eq. (16), 



C(t u t 2 ) = A'\tx 
x ln(|ti 



t 2 \ + 2B' -t 2 \ +a 2 )x 



i 2 |+a 2 )-a 2 ln(a 2 )] . 



(18) 



Evidently, this correlator is stationary and for \t\— 1 2 \ >> 
1, it has a power-law (with logarithmic correction) de- 
pendence on its argument as in Eq. (2) with H = 1/2. 
Thus for d — 2, the longitudinal motion x(t) of the par- 
ticle is a marginal fBm with H = 1/2. The analytical 
arguments leading to the persistence exponent 8 = 1 — H 
for fBm [15,16] can be easily modified to deal with the 
marginal case. This gives, P(t) ~ \Jh\{t)/t for large t, 
thus recovering the result in Eq. (1) for d = 2. 

In summary, we were able to show exactly that the 
longitudinal position x(t) of a particle in a (d + 1)- 
dimensional MdM model can be represented as a fBm 
with a Hurst exponent H(d) that varies with the dimen- 
sion d of the MdM model. Moreover, this exact con- 
nection allows us to use the known first-passage results 
for the fBm and derive analytically the asymptotic re- 
sults for the disorder averaged survival probability in the 
MdM model that were known before only via heuristic ar- 
guments and numerical simulations [6]. This method of 
connecting to fBm and the subsequent usage of the first- 
passage results for the fBm was first used to calculate 
exactly the first-return exponent of a fluctuating Gaus- 
sian interface to a fixed stationary configuration [16], a 
result which was subsequently verified in experiments on 
step edge fluctuations in Al/Si(lll) surface [17]. In this 
paper, we have shown that the same connection is use- 
ful in calculating the first-passage properties in a model 
with quenched disorder. Given that this technique is use- 
ful in two rather different problems, one expects that this 
method perhaps has much wider future applications. 

This paper has dealt with the persistence of a single 
particle in the MdM velocity field. It would be interesting 



to extend this study to a system of interacting particles. 
A simple physical example of such a system is a single 
polymer chain within the Rouse model where the beads 
on the chain are connected by harmonic springs [18]. The 
transport properties of a Rouse chain in presence of an 
MdM velocity field have been well studied in the liter- 
ature [19-21]. The techniques presented in this paper 
may be extended to study the persistence properties in a 
Rouse chain advected by an MdM velociy field [22]. 

It is a pleasure to thank S. Redner, P.L. Krapivsky and 
A. Comtet for useful conversations. 
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